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ABSTRACT 

It has recently been argued that single-collapse high-redshift models for elliptical 
galaxy formation can be rejected because they predict large numbers of very red 
galaxies at intermediate redshifts which are not seen in deep optical-infrared surveys. 
We argue, however, that this conclusion is premature since, while much effort has 
been invested in refining the predictions of hierarchical CDM models, only very sim- 
plistic models have been used to study the evolution of galaxies in other cosmogonies 
(e.g. isocurvature models). We demonstrate that the use of a more realistic multi-zone 
chemo-dynamical single-collapse model, yields colours at intermediate redshifts which 
are much bluer than inferred from the one-zone model, and indeed arc comparable 
to those predicted by hierarchical merging despite still allowing > 90% of the final 
stellar mass of elliptical galaxies to be formed in the first Gyr of their evolution. We, 
therefore, conclude that the one-zone model should be avoided to predict the colours 
of high-redshift galaxies and that the use of realistic multi-zone models allows the 
existence of ellipticals at high redshift, being their dismissal premature. 



1 INTRODUCTION 



The epoch of formation of elliptical galaxies remains a fun- 
damental and controversial issue in cosmology. While ob- 
servations show that at least some galaxies were already in 
place at redshifts z ~ 5 (Dunlop et al. 1996; Franx et al. 
1997; Dey et al. 1998; Hu et al. 1998), other recent data (Zepf 
1997; Abraham et al. 1998; Glazebrook et al. 1998; Kauff- 
man & Chariot 1998) have been interpreted as indicative of 
a general formation epoch at z < 3. In hierarchical models 
of galaxy formation dominated by cold dark matter (CDM), 
elliptical galaxies arise from the merging at low redshift of 
intermediate mass disks (e.g. Kauffmann, Chariot & White 
1996; Baugh, Cole & Frenk 1996). In contrast, in isocurva- 
ture CDM (Peebles 1998) (ICDM), or warm dark matter 
(Colombi, Dodelson & Widrow 1995) (WDM) models, el- 
lipticals are assembled at z > 3 in a major single collapse 
event. Zepf (1997) and Barger et al. 1998 have recently ar- 
gued that such models can be rejected because they predict 
large numbers of very red galaxies at intermediate redshifts 
which are not seen in deep optical-infrared surveys. However, 
this conclusion maybe premature since, while detailed pre- 
dictions (e.g. Baugh, Cole & Frenk 1996) have been made for 
hierarchical CDM models, only very simplistic models have 
been used to study the evolution of galaxies in ICDM or 
WDM models. In particular the physically unrealistic 'one- 
zone' model (Baade 1957; Arimoto & Yoshii 1987; Matteucci 
& Tornambe 1987) (or monolithic model), in which elliptical 
galaxies are formed in a single 10 8 -yr starburst predicts very 



red colours at all intermediate redshifts (Zepf 1997; Barger 
et al. 1998) . The aim of this paper is to explore a more real- 
istic multi-zone chemo-dynamical single-collapse model for 
predicting the colours of high-redshift elliptical galaxies. 



2 THE MODEL 

We have used the chemo-dynamical code developed by 
Friaga & Terlevich (1998) to compute a single collapse multi- 
zone model in which a single massive dark halo hosts bary- 
onic gas that will fall into the center of the dark potential 
and will subsequently form stars. The code assumes spher- 
ical symmetry and thus uses a ID hydro-dynamical solver 
that calculates the hydro-dynamical evolution and radial be- 
haviour of the ISM and stars during the early collapse and 
evolution of the elliptical galaxy. The galaxy is divided into 
100 zones from the core up to the tidal radius. We assume 
that the specific star formation rate (SF) follows a power- 
law function of gas density (p); SF(r,£) = vo[pl Po) 1 ^ 2 (here 
fo = 10 Gyr -1 , po is the initial average gas density inside 
the core radius of the dark halo, and r and t axe the spatial 
and time coordinates respectively). We have also included 
inhibition of star formation by inefficient gas cooling and 
expanding gas. A characteristic feature of these models is 
that several episodes of inflow and outflow occur simulta- 
neously at different radii (Friaga & Terlevich 1998). The 
detailed behaviour of these models of course depends on the 
precise parameters chosen. Although we will show here re- 
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is formed quickly enough to explain the strong metal lines 
observed in z ~ 5 QSOs (assuming that a young elliptical 
core hosts the QSO event), and the ICM is enriched rapidly 
enough to account for the lack of evolution of the ICM Fe 
abundance in rich galaxy clusters (Mushotzky & Loewen- 
stein 1997) up to z ~ 0.3. In addition, our multi-zone model 
does not display the extremely high luminosities (Lboi up 
to 10 15 L ) of the one-zone models, which have never been 
observed (Friaga & Terlevich 1998). Furthermore the model 
reproduces the observed spectral ages of red elliptical galax- 
ies found at intermediate redshifts (Dunlop et al. 1998). 
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Figure 1. The top panel shows the SF rate inside 10 kpc for our 
multi-zone model for three different masses. The primary epoch of 
star formation lasts for about 1 Gyr in each case, but, particularly 
in the most massive galaxy, low-level star-formation continues 
until the present epoch. Contrary to what is often claimed in 
the literature, a moderate extended period of star formation is 
consistent with the observed [a/Fe] ratios in elliptical galaxies. 
Since most of the stellar population is formed in a short burst 
of duration ~ 0.5 Gyr, the ISM is not contaminated by SNIa. 
The bottom panel illustrates this for two different masses, where 
the time evolution (not age) since the formation of the elliptical 
(compare with top panel) is shown. Note that 90% of the stars 
are formed at times smaller than 2 Gyr, thus the galaxy is [a/Fe] 
enhanced (about 0.4). Only stars born (less than 10%) after 2-3 
Gyr have moderate [a/Fe] ratios, but these do not contribute to 
the overall galaxy abundance. 

suits for specific cases, we believe that the critical feature 
is robust: the global star formation history in the galaxy 
naturally extends over a period of about 1 Gyr. 

Using the resulting distribution of the gas flow, the 
equations of chemical evolution are then solved. Once the 
structure of the ISM and the chemical abundances in each 
of the 100 zones have been computed, we use the spectral- 
synthesis code developed by Jimenez et al. (1998; 1999) to 
compute the spectral energy distributions of the multi-zone 
model at different radii and ages. Our models constitute a re- 
alistic family which successfully account for several observa- 
tional properties of elliptical galaxies (Jimenez et al. 1999, in 
preparation): (1) the observed metallicity gradients are nat- 
urally reproduced — the core has super-solar metallicity and 
the outer regions have sub-solar metallicities; (2) the [a/Fe] 
ratio is supersolar in the core of the galaxy, and also when 
averaged over the whole galaxy; (3) the observed positive 
mass-metallicity relation is reproduced; (4) galactic winds 
from elliptical galaxies can explain the amount of iron in 
the intra-cluster medium (ICM) of rich clusters of galaxies. 
The time scales for chemical enrichment are also consistent 
with the observations. For example, the high-metallicity core 



3 RESULTS 

One of the firm predictions of our model is that star for- 
mation in elliptical galaxies lasts longer than 10 s years. In 
the top panel of Fig. 1 we show the star formation rate in 
Mq yr _1 for three different galaxy masses. In general our 
elliptical galaxy model exhibits three star-formation stages: 
stage (1) during which the bulk of the stellar population is 
formed (half of the present day stellar population inside 10 
kpc is formed in 0.42, 0.40 and 0.30 Gyr for the galaxies 
with 10 10 , 10 11 , and 5 x 10 12 M , respectively); the dura- 
tion of this stage is well traced by the early, major peak, 
in the evolution of the star formation rate (see Fig. 1 top 
panel); stage (2), which corresponds to a extended period of 
star formation in the core (2-3 Gyr) as a result of the inner 
cooling flow. This episode has more modest star-formation 
rates than stage (1) and explains the blue colours of our 
model and thus the absence of passively evolving elliptical 
galaxies in deep optical and near-infrared surveys as claimed 
by Zepf (1997) and . Finally, after stage (2) is finished, our 
model evolves passively. However, very massive galaxies re- 
tain the inner (like e.g. cD galaxies in the center of clusters) 
cooling flow until the present epoch, which results in the 
occurrence of subsequent low-level bursts of star formation. 
For the intermediate mass galaxy the SF peaks at about 100 
Mq yr _1 . On the other hand, the most massive galaxies ex- 
hibit much higher maximum star-formation rates, typically 
several hundred MQyr -1 . Such a maximum star-formation 
rate is an order of magnitude smaller than that predicted 
by the one- zone model (~ 10000 Moyr -1 ), and is compara- 
ble to the maximum SF rates inferred for massive galaxies 
at 2 ~ 4 from sub-millimetre observations (Dunlop et al. 
1994). 

A recognized strength of the one-zone model is that a 
relatively short-lived burst of star formation seems to be 
required to avoid substantial contamination of the stellar 
population by Fe ejected into the ISM by supernovae SNIa 
(most of which explode after several 10 8 years), which makes 
it hard to reproduce the observed super-solar [a/Fe] abun- 
dances (where the a elements are: O, Ne, Mg, Si, S, Ar and 
Ca) in elliptical galaxies. Conversely, a problem in the hi- 
erarchical picture is that, due to the extended period over 
which significant star formation occurs, the resulting mod- 
els of elliptical galaxies struggle to exhibit super-solar [a/Fe] 
ratios, but rather tend towards the solar ratio due to SNIa 
contamination of the ISM. Since our multi-zone model has 
a somewhat more extended period of star formation (Fig. 1 
top panel) than the one-zone model, it might be expected 
that contamination by SNIa would also prevent the multi- 
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Figure 2. The time-evolution of V — K for three different masses and two cosmologies. We have plotted the time-evolution for the mass 
fraction contained inside a radius of 10 kpc. We have also assumed two formation redshifts: 3 (left line) and 10 (right line) for each mass. 
It is worth noting that in none of the cases the predicted V — K colours for the mass inside 10 kpc (the most realistic case, since at high 
redshift most of the observed light is inside 10 kpc), arc redder than V — K 7, two magnitudes bluer! than the expectations from the 
monolithic model as used in Zepf (1997). This strongly suggest that the dismissal of high redshift elliptical galaxies is premature. 



zone models from reproducing the observed [a/Fe] ratios. In 
fact this is not the case. In Fig. 1 (bottom panel) we show 
the time-evolution of [a/Fe] for two different masses and two 
radii. As expected, the predicted [a/Fe] ratios are initially 
very high due to the fact that SN II's dominate the chem- 
ical enrichment in the early evolution of the galaxy. Then, 
as later SNI explosions enrich the ISM with Fe, the aver- 
age [a/Fe] of formed stars decreases, but the star formation 
is actually switched off early enough to fix the [a/Fe] at 
super-solar values and it remains so until the present epoch. 
The reason for this is that the SNIa rate peaks during the 
late stage of continuing star formation (the maximum of the 
SNIa rate occurs at 0.85, 0.75, and 1.8 Gyr for the models 
with 10 10 , 10 11 , and 5 x 10 12 M Q , respectively), and by that 
time most of the stellar mass of the galaxy has been formed 
(> 90% — 95%). Therefore, the resulting galaxies still exhibit 
super-solar [a/Fe] ratios. 

The consequence of the relatively extended SF his- 
tory in this model is that elliptical galaxies do not be- 
come extremely red. Fig. 2 shows the redshift evolution of 
V — K for three different galaxy masses and two alterna- 
tive cosmological scenarios: an Einstein-de Sitter Universe 
(Q-m = 1.0, £7a = 0.0) and a vacuum-energy dominated Uni- 
verse (H m = 0.3, Qa = 0.7). These two example cosmologies 
have been chosen because they represent extreme cases of a 
young and an old Universe. In both cases we have adopted 
a value for the Hubble constant of 65 km s _1 Mpc -1 . In 
what follows we assume that all dark haloes collapse at a 
given redshift and then investigate whether, as the resulting 



galaxies evolve, they ever display observed optical-infrared 
colours redder than the limit considered by Zepf (1997) for 
the monolithic model, which translates into V — K < 8. For 
simplicity two formation redshifts are considered: z = 10 and 
z — 3. Note that at z > 3 the observed colours of galaxies de- 
termined by aperture photometry are usually roughly equiv- 
alent to their integrated colours, i.e. including the whole vis- 
ible radius. 

Fig. 2 shows that, even assuming a formation redshift of 
z — 10, the multi-zone model never produces colours as red 
as V — K — 8. Indeed it is only for the most massive galaxy 
model that a colour as red as V — K — 7 is ever achieved, 
and that requires the assistance of a very high formation 
redshift in an old (in this example A-dominated) universe. 
Thus, integrated within a radius of 10 kpc, the maximum 

V — K colour predicted by our multi-zone models is typically 

V — K ~ 6.5, comparable with the observed colours of the 
reddest known ellipticals at z > 1. Of course the maximum 
predicted colour is still bluer if a lower redshift of formation 
is assumed, and if z/ = 3, the maximum predicted colour 
peaks at V — K ~ 6 between z = 1 and z = 2. It is thus 
clear that even if all elliptical galaxies were to start form- 
ing stars at z — 10, they would never become as red as the 
one-zone model would predict. There is one and only rea- 
son for this: the extended period of star formation in the 
multi-zone model. Note that these blue colours are indepen- 
dent of the suite of spectrophotometric models chosen to 
compute them, since the effect of choosing different codes 
will be at the most 0.3 magnitudes (see e.g. Jimenez et al. 
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Figure 3. The abundance of faint galaxies in the field, collected 
by Zepf(1997), is compared with the predictions of the multi-zone 
model for formation redshifts of 10 (dashed line) and 3 (continu- 
ous line). For a colour cut-off of V — K > 8 the multi-zone model 
predicts zero galaxies (even for zj = 10) in superb agreement with 
the observations. Moreover, even at a much lower colour cut-off 
(V — K > 6), the multi-zone model predicts a very low surface 
density of red objects, in agreement with the observations (now 
displayed as a shadow region). 



1999), totally irrelevant for the present discussion. In other 
words, small differences in rest frame colours have no effect 
in observed colours. A direct comparison with the observed 
data collected by Zepf (1997) on the abundance of faint red 
galaxies is shown in Fig. 3. The top two panels show that 
the multi-zone model predicts zero galaxies with V — K > 8, 
in superb agreement with the observations. An even more 
robust prediction of the multi-zone model is that the space 
density of galaxies with V — K > 6 (two magnitudes below 
the colour cut-off chosen by Zepf(1997)) will be smaller than 
five galaxies per arcmin 2 again in superb agreement with the 
data. 

Finally, we consider the star-formation rates estimated 
from the single collapse multi-zone model in the context of 
the cosmological evolution of global star formation rate in 
the Universe. Several authors have recently attempted to 
determine the global level of star formation in the Uni- 
verse at redshifts higher than two (Pei & Fall 1995; Dun- 
lop 1998; Pettini et al. 1998; Hughes et al. 1998); the main 
uncertainty in this case is the precise impact of dust at 
high redshift. To predict the star-formation density con- 
tributed at high redshift by young ellipticals from the sin- 
gle collapse multi-zone model, we have assumed that the 
number density of massive ellipticals at high-redshift is 
the same as at the present-day (Terlevich & Boyle 1993), 
and that they all collapse at the same redshift. There- 
fore we adopt a Schechter luminosity function given by 
$>{M B )dM B = $*10°- 4(m b- Mb) ° exp[-10 a4(M S- M B>] with 
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Figure 4. The figure displays the global SF in the Universe as- 
suming that all elliptical galaxies form at a certain redshift: 4.5 
(solid line) and 6 (dashed line). Also plotted are several observa- 
tional indicators of the cosmological evolution of star-formation 
density, derived from Lyman-limit galaxies in the HDF(asterisks, 
the point at z = 4 has been also plotted (highest point) includ- 
ing the correction from incompleteness due to colour selection 
in the HDF (C. Steidel, private comm.)), from recent deep sub- 
millimetre surveys (triangle), from the cosmological evolution of 
luminosity density from radio-sources (dotted line), and from the 
chemical evolution of metals in the Universe (shaded region). It 
transpires from the plot that, even if all elliptical galaxies were 
formed at redshift 6, their predicted contribution to global star- 
formation density would still be consistent with the existing data 
and the other external indicators and there would even be room 
for the bulges of spiral galaxies to be formed at high redshift. 



3 h 3 Mpc~ s with 
—22.5 adopted 



a = 0.23, M B = -21 and <F = 2.9 x 10" 
the same power law extension for Mb < 
in Friaga & Terlevich (1998). Integrating this luminosity 
function for galaxies brighter than 0.1L„, it is possible to 
predict the global star-formation density of the Universe 
at any redshift using our multi-zone models for different 
galaxy masses. Fig. 4 shows, for an Einstein-de Sitter Cos- 
mology, and two different formation redshifts (z/ = 4.5; solid 
curve, and Zf = 6; dashed curve), the predicted contribu- 
tion of young ellipticals to the global star-formation density 
at high-redshift. For comparison, we also show several ob- 
servational indicators of the cosmological evolution of star- 
formation density, derived from Lyman-limit galaxies in the 
HDF (Pettini et al. 1998), from recent deep sub-millimetre 
surveys (Hughes et al. 1998), from the cosmological evolu- 
tion of luminosity density from radio-sources (Dunlop 1998), 
and from the chemical evolution of metals in the Universe 
(Pei & Fall 1995). Fig. 4 demonstrates that our models do 
not conflict with the existing data, and other external indica- 
tors of star-formation density at high redshift. Even galaxy 
formation redshifts as high as 6 appear perfectly consistent 
with the available data. Note that we have used an extreme 
case in which all galaxies form at the same redshift. If ellipti- 
cal galaxies form over a redshift range between 3 and 10 with 
a small fraction of them being formed at redshift higher than 
6, then there would be even room to form bulges of spirals 
at a high redshift. 
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4 CONCLUSIONS 

We conclude by re-emphasizing that the hierarchical picture 
and the multi-zone model considered here both predict a 
complete absence of objects with colours as red as V — K ~ 7 
in deep optical and near-infrared surveys (unless, of course, 
the intrinsic colours of the galaxy in question are heavily 
distorted by dust). Furthermore the multi-zone model does 
display some late episodes of low-level star formation that 
could account for the observed properties of some ellipti- 
cal galaxies (Abraham et al. 1998; Glazebrook et al. 1998) 
in the HDF which appear to show evidence of some recent 
star-forming activity and it does not seem implausible that 
field and cluster ellipticals may follow a different evolution- 
ary path. However, in the multi-zone model considered here, 
> 90% of the final stellar population of the final elliptical 
galaxy is formed during the first Gyr of its evolution. This 
provides a natural explanation for the homogeneous prop- 
erties displayed by present-day elliptical galaxies, without 
violating the colour limits derived by Zepf (1997) at inter- 
mediate redshifts. The lack of very red galaxies at high red- 
shifts is therefore not a strong test of hierarchical structure 
formation. Better probes of such theories come either from 
the homogeneity of the colour-magnitude relation (Bower, 
Kodama & Terlevich 1998) , or from a direct counting of lu- 
minous galaxies at high redshifts. At present, there is some 
evidence that the galaxy population was not created in its 
entirety at very high redshifts (Kauffmann & Chariot 1998) , 
but single-collapse models for the spheroidal population can 
not yet be ruled out (e.g. Totani & Yoshii 1998). Larger com- 
plete samples with HST morphologies will be needed before 
we can claim to have observed merging at work in modifying 
the galaxy population. 
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